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Abstract. Using supersymmetry techniques analytical expressions for the average of 
the fidelity amplitude fe{T) — {il;{0)\exp{2TnH^T) exp{—2TnHoT)\tp{0)) are obtained, 
where He — Hn + {^/e/2Tr)V, and Hq and are taken from the Gaussian unitary 
ensemble (GUE) or the Gaussian orthogonal ensemble (GOE), respectively. As long 
as the perturbation strength is small compared to the mean level spacing, a Gaussian 
decay of the fidelity amplitude is observed, whereas for stronger perturbations a change 
to a single-exponential decay takes place, in accordance with results from literature. 
Close to the Heisenberg time t = 1, however, a partial revival of the fidelity is found, 
which hitherto remained unnoticed. Random matrix simulations have been performed 
for the three Gaussian ensembles. For the case of the GOE and the GUE they are in 
perfect agreement with the analytical results. 
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1. Introduction 

The concept of fidelity has been developed as a tool to characterize the stability of a 
quantum-mechanical system against perturbations pp. Originally fidelity was introduced 
as the squared modulus of the overlap integral of a wave packet with itself after the 
development forth and back under the influence of two slightly different Hamiltonians. 
Let Hq be the unperturbed Hamiltonian and 

He = Ho + ^V (1) 

the perturbed one. This somewhat unusual definition of the perturbation strength e has 
been applied for later convenience. Then the fidelity is given by 

F,(r) = |(V'(0)|exp(27rziJ,r)exp(-27r2ifor)|^(0))|' , (2) 

where ip{0) is the wave function at the beginning, often chosen as a Gaussian wave 
packet with minimum uncertainty. It is assumed that Hq has mean level spacing of one, 
and thus r is given in units of the Heisenberg time. The variance of the off-diagonal 
elements of V is chosen to be one. In all what follows it is assumed that e is of the order 
of one thus guaranteeing that the shift of the levels due to the parameter variation is of 
the order of the mean level spacing. 

Depending on the strength of the perturbation one can discriminate roughly three 
regimes. In the perturbative regime, where the strength of the perturbation is small 
compared to the mean level spacing, the decay of the fidelity is Gaussian. As soon 
as the strength of the perturbation becomes of the order of the mean level spacing, a 
cross-over to exponential decay is observed, with a decay constant obtained from Fermi's 
golden rule j21 E] • For very strong perturbations the decay becomes independent of the 
strength of the perturbation. Here the decay is still exponential, but now the decay 
constant is given by the classical Lyapunov exponent ji]. It has been proposed by 
Pastawski et al ^ to look for such a behaviour in a spin-echo experiment on isolated 
spins coupled weakly to a bath of surrounding spins [H]. 

A paper of Gorin et al jZj is of particular relevance for the present work. The 
authors calculated the Gaussian average of the fidelity amplitude in the regime of small 
perturbations using the linear-response approximation, 

/,(r)~l-6C(r). (3) 

where C(r) is given by 

r2 r '■^ 



(3 2 



JO 



CW = ^ + o- / / MOdt'dt, (4) 



and 1 — b2,j3{T) is the spectral form factor. (3 is the universality index, i. e. = 1 for the 
Gaussian orthogonal ensemble (GOE), (3 = 2 for the Gaussian unitary ensemble (GUE), 
and (3 = 4: for the Gaussian symplectic ensemble (GSE). For an explicit calculation 
knowledge of the spectral form factor is thus sufficient. By an exponentiation of the 
above formula, 

/,(r)~e-^M, (5) 
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the authors were able to describe the cross-over from Gaussian to exponential decay with 
increasing perturbation strength quantitatively. The Lyapunov regime is non-universal 
and thus not accessible in a random matrix model. 

It is obvious that the linear-response approximation must break down for large 
perturbations. In the present work supersymmetry techniques are applied to calculate 
the ensemble average of the fidelity decay. Since this calculation is non-perturbative, 
the results hold for arbitrary values of the perturbation strength. We shall see that 
the calculation reveals an important generic feature, which is unaccessible by any 
perturbative approach. 



2. The model 



In the present paper we shall discuss the ensemble average of the fidelity amplitude, 
since for this quantity the calculation is much easier than for the originally introduced 
quantity Q. Since both the unperturbed Hamiltonian and the perturbation are taken 
from the Gaussian ensembles, the choice of the initial wave packet 'ip{0) is irrelevant. 
The ensemble average may thus be written as 

/,(r) = l(Tr[e2-^-e-2-^-]>, (6) 

where it is assumed that the Hamiltonian has been truncated to a finite rank A^. 

The Hamiltonian introduced in equation has the disadvantage that the mean 
density of states changes with e. The more it is somewhat inconvenient for the present 
calculation that the variances of the matrix elements of Hq and V differ. We therefore 
adopt a slightly different parameter variation, 

= Hq cos (f) + Hi sin (j) . (7) 

It is assumed that the matrix elements of Hq and Hi have the same variance, have zero 
average, {{Hq)^) = {{Hi)ij) = 0, are uncorrected, ((ifo)ii(-f^i)A:/) = 0, and are Gaussian 
distributed, 

p({(i7„),J)~exp(^-^Tri72^ , n = 0, 1 , (8) 

where 

From random matrix theory it is known that the corresponding ensemble averaged 
density of states is given by Wigner's semi-circle law with a value of one in the centre 
of the circle. 
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(see e.g. reference |Bj). The resulting variances of the matrix elements of Hq and Hi 
are given by 

It follows for the GUE 

m^UH^Aki) = -cos{<P2- 'Pi) . (12) 

A similar expression is obtained for the GOE. Since only the difference of 02 and 0i 
enters expression (fT^ . we may assume without loss of generality 0i = —02 = 0/2. 
Ansatz ((7j) for the parameter variation is obtained from equation ^ by means of the 
substitutions 

TT 

iJo — ^ COS iJo , ^ ^ ' ^^^^ 

where 



tan0=^^. (14) 

is thus of 0{^^), and e is given in the limit of large A^ by 

e = 4Ar02 . (15) 

Only terms up to 0{(f?) will survive the limit A^ — cxd as we shall see later. The details 
of the parameter dependence are irrelevant. 

With all these substitutions equation (jH)) may be transformed into 

/,(r) = j dEi dE^e^-^''^-'^-^^ {E^, E^) (16) 

where 

fi.(B,B,).^(T.(^— -i^)). (17) 

with E± = E ±zri, and the abbreviations c = cos(0/2), s = sin(0/2). Using standard 
supersymmetry techniques 0, this can be written as 

-Re (El, E2) 

^ Jl\ I ^[^] XI (^"^™ ~ CU){y*myn - VmVn) 

^ n,m 

^ ^-ix^Ei-cHo~sHi)yL ^'i.yHE2-cHo+sHi)y 

= ^/ d[x]d[y] $^(x:x„-a«.)(l/;nZ/n-r/;;r^n)e-4''^^--^^^^^] 

n,m 

^ ^gjc[xtHox-yt/foy]^ ^gis[xt//ix+yt/fiy]^ ^^^-^ 

where x = (xi, ^1, . . . , xn, ^jv)^, y = (yi, • • • , l/jv, Vn)'^ , and 

= n ' ^ [^] = n ^y*^ ^^'^ ■ ( 1 ^) 
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We adopt the usual convention and use latin letters for commuting, and greek ones 
for anticommuting variables, respectively. Equation (fTHj) is still true for all Gaussian 
ensembles, but now we have to discriminate between the GUE and the GOE. 



3. The GUE case 



Using equation (jH)), the calculation of the Gaussian average over Hq and Hi is elementary. 
The result for Hq may be expressed as 



g*c(xt//ox-ytHoy)^ = e"^ 



(20) 



5 = E 



V^n' Sn' Vn) ' 



(21) 



where 

-Vn 

V / 

Whenever supermatrices are involved, traces and determinants are to be interpreted 
as super traces and determinants, respectively, in the definition of reference [H] . In short 
hand notation equation may be written as 

5 = L^z„zt, (22) 



where 

Vn 
\Vn J 

and I2 is the two-dimensional unit matrix. Introducing the notation 

Saa Sar 
Sra Srr 



(23) 



S 



(24) 



where each Sij is a 2 x 2 matrix, and the indices 'A', 'R' refer to the 'advanced' and 
'retarded' components, respectively, the sum entering equation (fTH|l may concisely be 
written as 



^{x^Xn, - Cnim){y*myn " VmVn) = -Tt{SaR CT SbA Ct) , 



where 



n,m 



a 



Next, a Hubbard- Stratonovich transformation is applied to equation ()20|). 



(xtHox-yt/foyl 



(i[M]e 



(25) 



(26) 



(27) 
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where U is the supermatrix 



Ura Urr 

with the 2x2 components 



I I , 1,3 = AR. (29) 



For the integrals in equation (j^7j) to be well-defined, the Uij integrations have to 
be performed from — oo to oo, whereas the Uij integrations are from —too to too. (In 
literature usually an additional factor of t is introduced in the lower right corner of the 
matrix ()29|1 to avoid integrations along the imaginary axis.) 

In the same way we obtain 

^^.(xtH.x+yt^,,)^ = I d[v]e~^^'''^''^''\ (30) 

where 



and 



with 



V.j= "^'^ , ^,J = A,R. (33) 

V J 

Collecting the results we obtain from equation (jl8|) 

d[u] d[v]e ^^f< Tr a a 



cm J ' ' ' ' \dUAR OUra 

Now the X, y integrations can be performed resulting in 



where 

El I2 - cUaa - sVaa cUar - sVar 

-cUra - sVra -E2 12 + cUrr - sVrr 



\M\ 



(35) 
(36) 
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Introducing the notation E1/2 = E ± E/2, this can be written as 



\M\ 



EU-CU+ (^-sV ]L 



(37) 



where I4 is the 4x4 unit matrix, and L has been given in equation (j'23j) . Substituting 



U 



U+{j-sV\L 



(38) 



the V integrations can be performed yielding 

Re (El, E2) 

ih T^[UlA + U%B+mAA-URR)+^{UARU[iA) 



- I d[u]e 



X Tr 



aT:^^(r ] lEU 



OUar OUra 



-N 



NX} COS2 : 



d[u] Tr {Uar cr Ura(t) 



^ ^-2t^'^[mAA-Unn)+2{^-l)iUARURA)]^-Trg(U) 



where 



9iU) 



El 

2\f 



+ Arin(El4 - U) 



(39) 



(40) 



The second equation (jH^ has been obtained by an integration by parts. 

Equation (jH^ is still exact, but now the limit — cxd is performed. Since A/3 is of 
0{N) (see equation 0), U and E in equation (jlOj) . too, must be of 0{N). E, on the 
other hand, is of the order of the mean level spacing and thus of 0{1). Furthermore, 
is of 0{^^). Consequently Tr (?([/) is of 0{N), whereas all other terms entering the 
integral (j^ are of 0(1). 

This suggests to diagonalize U, 

U = T-^UdT, (41) 

and perform the integrations over the elements of the diagonal matrix Ud by means 
of the saddle point technique. The saddle points are obtained from the zeros of g'{u), 
whence follows 



N / nE 



ua/r 



TT 




(42) 



D)s 



The plus and the minus sign belong to the advanced saddle point ua, and the retarded 
one ur, respectively. The matrix Ud at the saddle point is thus given by 

UAI2 
URI2 

The matrix T diagonalizing U may be parameterized as 

VI + 5C iB 
-iC Vl + Cfi 



T 



(43) 



(44) 
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where B and C are 2x2 supermatrices. Inserting equations (j43j) and (|44p into equation 
|^T|) . we obtain for the matrix U at the saddle point 



^ , uaU + i^BC ABVTTCB 



where 



2N / /7r^\^ 2A^ 

In the last equation we used expression (|TO|l for the mean density of states p. 
We are now left with 

Re (El, E2) 

1- /Tr (f/^« a UnAcr) e"^ ^[^(^--^-)+*^(^-^-)] \ 



iVA^X /' 

(47) 

where only terms in surviving the N ^ 00 limit have been taken. The brackets 
denote the average over the angular variables entering the matrix T, taken at the saddle 
point. Using equation we obtain for the quantities entering on the right hand side 
of equation (P7|) 

Tr {Uaa - Unn) = 2tA Tr BC , (48) 
Tr UarUra = Tr [BC + (BC)^] , (49) 
TiiUARcrURAcr) = Tr(CVl + BCaBVl + CBa) . (50) 
The matrices B and C are diagonalized by means of the transformation 

B = PBdQ-^ , C = QCdP-^ , (51) 



where 



Bo^i' Co^l': ^ I. (52) 




a* v/TT^ ; ' ^ \ (3* virm ) ' ^ ^ 

(see e.g. Chapter 10 of reference |10i). Inserting these expressions into equations (jlSjl 
to (jHUjl . we obtain, 

Tr {Uaa ~ Urr) = 2mpE {tf + W) , (54) 

^ Tr UarUra = ^p' [tt* + tt* + [trf - [Wf] , (55) 
i:i{UARaURAa) = A^ Tr(Ci3 x/l + BdCdctpBdsJI + CdBdctq) , 

(56) 
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where expression (P) for \p and expression (fT^ for e were used, and 

^ y -2a* -l-2a*a / ' ^ ^ 

n-i n ( 1 + 2(3 \ 

= ^ = ( -2/3* -1 - 2/5*/3 j ■ ^^'^ 

For the calculation of the average (j47j) over the angular variables the 'surface volume' 
element is needed, 

dt dt* dt di* da da* dp dp* . . 

"""l = (iFTWf 

(see again reference JO])- The integral over the anticommuting variables is easily 
performed. Only ap and aq depend on the variables a, a*, and P,P*, respectively, 
and the corresponding integrals reduce to 



da da* ap ~ I2 , f dp dp* ctq ~ I2 , (60) 
whence follows 

jdada- dl3dl3'MU,,.U^.) ~ Tr(U,nU,,) 

~ A^Tr [BC^{BCf] 
~ [tt* +tt* + {tt*f - {tt*f] . (61) 
Collecting the results, we obtain from equation (jTTj) 

1 f^\" f . ,*tt*+tt* + {tt*f-{tt*f 



R, {El, E2) ~ ^ (^Y^ j dt dt* dt dt 



{tt* +tt*f 

The t, t* integration is over the whole plane, whereas the t, i* integration is restricted 
to the unit circle tt* < 1. Introducing polar variables, we obtain 

00 1 

U.iE.,,,.^jd.jdy^±l±^ 



Inserting this result into equation (fTI)jl . and introducing E = {Ei + E2)/2 and 
E = El — E2 a.s new integration variables, we get, fixing the constant of proportionality 
by the condition /e(0) = 1, 

00 1 

f.{r) = l fdEp' [dx [dy^ + ''-y 



N J J J X + y 





X S[t - p{x + y)]e-^p'(-+y)(^+--y) . (64) 




Figure 1. Fidelity amplitude /e(T) for the GUE (solid lines) and the GOE (dashed 
lines) for different values of the perturbation strength e. r is given in units of the 
Heisenberg time. 



The E integration is nothing but an energy average. Restricting the discussion to the 
band centre, we may discard this average and obtain 



1 



r 



Min{T,l) 



/,(r) = - / dyil + r- 2y)e-t-(i+-2,) _ (gg^ 



The integral is easily performed with the result 

s(f r2) - rs'(f r2 
^ e 2' "[s(er) - is'(fr)] , r>l 

where 



sinh(x) 

s{x) = , 67 

X 

and s'{x) denotes its derivative. 

Equation (jUBj) is the central result of this section. It gives an analytic expression for 
the GUE average of the fidelity amplitude for arbitrary perturbation strengths, /^(t) 
and its first derivative are continuous, but the second derivative shows a discontinuity 
at r = 1. A similar situation is known for the spectral form factor, where, however, for 
the GUE already the first derivative is discontinuous. 

The solid lines in Figure ^ show the GUE fidelity amplitude for different values of 
the perturbation strength e. For e -C 1 the fidelity decay is predominantly Gaussian. 
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For e = 1 and small times r an exponential decay is found, with a cross-over to Gaussian 
behaviour at r = 1, both observations in accordance with results known from literature. 
For e ^ 1 the fidelity decay is exponential for short times. Close to r = 1, however, 
there is a conspicuous partial revival of the fidelity which had not been reported before, 
as it seems. For still longer times the decay becomes Gaussian again. 
In the limit of small perturbations equation reduces to 



Mr) 



l-f(^r + ^j, r<l 
l-f(i + r2), r>l 
This is in complete accordance with the results obtained by Gorin et al jTj. 



(6^ 



4. The GOE case 

The first steps in the calculation of the ensemble average of the fidelity amplitude for 
the GOE are the same as for the GUE. Equation for S remains correct, but now S 
is an 8 X 8 matrix with z„ given by 



(69) 



]_ c 

V2 Vn 
Vn 
Vn 
\V*n ) 

In taking the adjoint of z„ one has to consider that the complex conjugate of 
the complex conjugate of an antisymmetric variable is defined as (a*)* = —a, see the 
appendix of reference 

Up to equation the further procedure is nothing but a step-by-step repetition 
of the calculation for the GUE case. The main problem for the GOE case arises from 
the diagonalization of matrix ?7, 

U = T-^UdT , (70) 

see equation ()4H). Not all of the matrix elements of S are different, as is evident from 
its definition, with the consequence that 5* obeys a number of symmetries which are 
inherited by the matrix U. The matrices T have to be chosen such that all symmetries 
are conserved. It is a highly non-trivial task to find the best parameterization for the 
matrix elements of T obeying these constraints. Fortunately this problem has already 
been solved by Verbaarschot, Weidenmiiller, and Zirnbauer in their disseminating work 
[Oj. We just cite their essential results: 
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Equations still hold, but now Bd and Cd are equal and given by 

fh ■ ■ ■\ 

■ it 

■ ■ it 



B 



D 



Cd 



\ 



(71) 



The parameterization of the matrices P and Q is complicated, and is given in the 
appendices of reference p. For the present purpose it is sufficient to note that the 
angular averages over the matrices crp = PaP~^ and aq = QaQ~^ (see equations (|H7jl 
and()58|)) up to a constant factor again yield the unit matrix, as can be shown by explicit 
calculation. All formulas of section El can thus be applied directly to the GOE situation. 
The surface volume element for the only remaining variables ti,t2,i is given by 



d[n] 



-dti dt2 dt 



(72) 



t?)(l+ti)(t?+t2)2(t2 + p)2 

and the integrations are from to oo for ti and t2, and from to 1 for i. 

Collecting the results, and proceeding in exactly the same way as for the GOE case, 
we finally end up with 

oo oo 1 

y\ 



Mr) 



dx dy dz 



, 1 — z)z\x 



Substituting u 

Mr) 



^/x{l + x)y{l + y){x + z)'^{y + z)^ 

u u u 

X [x(l + x) + y{l + y) + 2z{l - z)] 
X + y)/2 and v = {x — y)/2, we obtain 

T U 

vdv (r — ■u)(l 



(73) 



du 

~1) 



T + u) 



Max(0,- 

X \{2u + l)r - + ^2]e-iI(2«+i)r-r2W] ^ (74) 

where the constant of proportionality again was fixed by the condition /^(O) = 1. For 
e = the right hand side of equation (f7H) must be one by construction, but it is not 
straightforward to show this explicitly. Since the corresponding calculation may be of 
some interest, it is reproduced in [Appendix A[ 

Equation ([71)) gives an explicit expression for the fidelity amplitude for the GOE 
case. It is not yet suited directly for a numerical integration, since the integrand contains 
a number of singularities. But it is not difficult to remove them by suitable substitutions 
of integration variables. This is done in [Appendix B 



The dashed lines in Figure ^ show the results of the calculation for the same e 
parameters as before. We notice that the partial recovery of the fidelity close to r = 1 
is still present for large e values, but is considerably less pronounced than for the GUE 



case. 
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5. Numerical simulations 




In this section we present random matrix simulations to affirm the analytical findings 
for the Gaussian orthogonal and unitary ensembles. Further we show numerical results 
for the Gaussian symplectic ensemble which has not been treated analytically. 

In our simulations the Hamiltonians Ho and Hi are random matrices of dimension 
N X N with variances of the diagonal and off-diagonal elements given by equation (|TT|). 
To calculate the fidelity amplitude, we write expression © as 

Mr) = 1 (Tr ^/^,e^-^^i^,-^/^oe-^-^°"^i^o^ 

'-^'-^'''W"^ , (75) 

where Hq = Rq^HqRq and H^ = R^^H^R^ are diagonal, and R = R'^^Rq. 

In the numerical simulations the trace in equation ()75|1 was restricted to 20 percent 
of the eigenvalues in the centre of the spectrum where the mean level density is still about 
constant. The average was taken over up to 8000 random matrices for Hq, and for each 
of them over 50 random matrices for Hi. For larger values of the perturbation strength 
e it became more and more important to choose the dimension of the matrices large 
enough to avoid finite-size effects. = 500 proved to be sufficient for e < 10. 

The results for the three Gaussian ensembles are shown in Figure El For the GOE 
and the GUE the numerical simulations are in perfect agreement with the analytical 
result for all e values shown. For comparison, the fidelity amplitudes in the linear 
response approximation [Z] (see equations (jHI) and (jl}) are shown as well. For small 
perturbation strengths and small values of r, the linear response result is a good 
approximation, but the limits of its validity are also clearly illustrated. In particular, it 
does not show any indication of the recovery near r = 1. 



6. Discussion 



This work extends the results by Gorin et al ^ to the regime of strong perturbations 
using supersymmetry techniques. An intuitive explanation for the surprising recovery 
of the fidelity amplitude at the Heisenberg time can be given in terms of the Brownian- 
motion model for the eigenvalues of random matrices. The behaviour of the fidelity 
amplitude has its direct analogue in the Debye- Waller factor of solid state physics (see 
reference It is stressed that our result is generic and not restricted to random 

matrix systems. For instance, the fidelity recovery has recently been observed in a 
spin-chain model by Pineda et al |12j . 

The results of the present work may be easily extended to all situations, where 
the Gaussian averages (see equation (fTK|l ) lead to expressions allowing a subsequent 
Hubbard-Stratonovich transformation. This is, e.g., the case, if Hq is taken from the 
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1 2 3 4 5 



T 




1 2 3 4 5 



T 




T 



Figure 2. Fidelity amplitude /^(t) for the GOE (a), the GUE (b), and the GSE (c) for 
e = 0.2, 1, 2, 4 and 10. The solid lines show the results of the numerical simulations, 
and the dotted lines those of the linear response approximation. For the GOE and 
the GUE the numerical results are in agreement with the analytical results within the 
limits of the line strength. 
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GOE, and Hi is from the GUE, or is purely imaginary antisymmetric. Perturbations, 
where the diagonal is zero, are of particular interest, since in such a situation the decay 
of fidelity freezes [13113. 

There might be still another application of the formulas derived in this paper. If 
/e(r) is expanded into a power series of e, the linear term can be expressed in terms of the 
spectral form factor, i. e. the Fourier transform of the two-point correlation function [7j. 
In a similar way the coefficients of the nth power of e depend on all fc-point correlation 
functions up to A; = n. fe{T) may thus be used as a generating function to obtain these 
terms in a simple way. 
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Appendix A. Proof of /o(t) = 1 



Substituting x = \/u^ — v'^, we obtain from equation (fTijl 

T U 

dx {t — u){l — T + u) 



fe{T) =2 j du 

Max(0,T-l) 



X [{2u + l)r -r'^ + u'- a;2]e-f . (A.l) 



Specializing to e = 0, and applying the substitution 

x={2u + 1) 
this may be written as 



2 V y2^TT 



where 



/o(r) = 2 j du{r-u){l-T + u) 

Max(0,T-l) 

X [{u" + 2™ + r - r')/2 - {2u + l)h] , (A.2) 



u 

dz ~^ ^ , (A.3) 

[(r2-M2)(2z+l) + (2M+l)z2]2' ^ ' 

u 

/2 
dz . . (A.4) 
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The same equations can be found already in reference ^H] • The latter two integrals 
can be performed elementary and yield 



— uq u(u + 1) 
h = ^ o o arctan - — —7^ h 



2p3g3 

1 



2pq^ 



arctan ■ 



{u + l)p 2r2p^g2 
uq u{u + 1) 



where 



P 



u + 2r2(2u + l)g 



V(« + l)2-r2 



2 ' 



(A.5) 
(A.6) 

(A.7) 



Let us denote the part of f^ij), depending in the integrand of the arctan terms, by 
fQ[j)- It may be transformed as 



du 



p2 q2 



(2m+ 1)- 



Max(0,r-1) 

X r — n l + n — r — arctan — 

^ ^pq {u + l)p 



du 



Max(0,r-1) 



+ 2tu + t -t"^ + 2u{t + 1) + r - + 1 
{t-u){t + u) (m + 1 - r)(M + 1 + r) 



X 



(r — + 1 — r) 

~ T"; ^ cLrctjcLll ^ — 

{t + u){u+1 + t) {u + 1)p 



du 



2 + 



2u + l 



1 



Max(0,T-l) 



X 



(r — + 1 — r) 
(r + + 1 + r) 



arctan 



u — r u + T — r u + l + r 



2 + (2m + 1) 



d_ 

du 



Max(0,T-l) 



(r — m) (m + 1 — r) 
{t + u){u + 1 + t) 



arctan ■ 



uq 



du 



r 



+ 1 



Max(0,r-1) 



{t + u){u + 1 + t) 



arctan ■ 



[u + l)p 

uq 
u + l)p 



+ 



du{2u+ l)^ 



Max(0,T-l) 



{t — u){u + 1 — t) d ^ uq 
(r + 1 + r) du {u + l)p ' 



(A. 



where in the last step an integration by parts was performed. The terms in the first 
row cancel, and only the term in the last row survives. Performing the differentiation. 



Recovery of the fidelity amplitude for the Gaussian ensembles 
we have 

l + 3M(M + l)-r2 



/o(r) 



du- 



{u + t){u + 1 + t) ■ 

Max(0,T-l) 

Collecting the results we are left with 
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(A.9) 



/o(r) 



du 



1 



(M + r)(M + 1 + r) 



[U^ + 2TU + T- T^)^-^ 



Max(0,r-1) 

, u(u + l)(r^ — M^) / X ' 

+ (2^ + 1) ^ (2^ + 1)^2 ^ +l + 3^^+l)-r 

r 

+ 3mV - + + 3Mr + u + t 



1 

r 



Max(0,T-l) 



{u + t){u + 1 + t) 



^ y du(2M - r + 1) 



Max(0,T-l) 
1 , 

— — r + 1) 
r 



Max(0,T-l) 



q. e. d. 



(A.IO) 



Appendix B. Transformation of the integral (|74p 

To turn equation (f7^ into a form being suited for a numerical calculation, we start with 
equation ()A.1|) by substituting m = r sin and obtain 

7r/2 

/e(r) = 2 y r cos (/) r(l — sin </))[! — r(l — sin (/))] 



Max(0,arcsin[(r-1)/T]) 
T sin (/) 

X 



dx 



a/o;^ + 2r sin + 1 (x^ + cos^ 
where 

Z = (2r sin + l)r — cos^ (j) — . 
Next we substitute x = xr cos cf) with the result 

7r/2 

m = 2 1 

Max(0,arcsin[(r-1)/T]) 



1 — t(1 — sin I 



tan f 



X 



1 + sin ( 
dx 



x^r^ cos^ + 2r sin + 1 (1 + x'^Y 



:Ze~^ 



(B.l) 



(B.2) 



(B.3) 



Recovery of the fidelity amplitude for the Gaussian ensembles 
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where 

Z = 2rsin0 + 1 - rcos20(l + £2) . (B.4) 
After the final substitution x = tan a we end with 



,(.) ^ 2 / 



7r/2 

, 1 — t(1 — sin I 



1 + sin ( 

M£ix(0,arcsin[(r-1)/T]) 



/ , COS a [(2r sin </) + 1) cos^ a — r cos^ i«i _s^7 ,^ , 
X / da =e 2 , (B.5) 

y cos^ sin a + (2r sin 0+1) cos^ a 



where now 



cos^ 



2rsin0 + l-r . (B.6) 



cos^ a 

The integrand of the double integral ()B.5|1 behaves well everywhere, and the 
numerical integration does not pose problems any longer. 
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